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ABSTRACT 

We present a study of 15 new brown dwarfs belonging to the ~ 7 Myr old 25 Orionis 
group and Orion OBla sub-association with spectral types between M6 and M9 and 
estimated masses between ~ 0.07 Mq and ~ 0.01 Mq. By comparing them through 
a Bayesian method with low mass stars (0.8 ^ M/Mq ^ 0.1) from previous works 
in the 25 Orionis group, we found statistically significant differences in the number 
fraction of classical T Tauri stars, weak T Tauri stars, class 11, evolved discs and 
purely photospheric emitters at both sides of the sub-stellar mass limit. Particularly 
we found a fraction of S.Otf.g % low mass stars classified as CTTS and class 11 or 
evolved discs, against a fraction of 33.3lg°g® % in the sub-stellar mass domain. Our 
results support the suggested scenario in which the dissipation of discs is less efficient 
for decreasing mass of the central object. 

Key words: stars: low-mass, brown dwarf, open clusters and associations: individual 
(25 Orionis) 


1 INTRODUCTION 

After two decades of studying the properties of young brown 
dwarfs (BD), there is an agreement that these objects and 
very low mass stars (VL MS) have similar formation pro¬ 
cesses (e.g. lLuhmarill2012h although some issues about their 
early evolution are still matter of intense research. For in¬ 
stance, it is well accepted that the number fractions of 
VLMS and BD harbouring primordial discs drop off dur¬ 
ing t he first ~ 10 Myr of their evolution (e.g. iLuhmai] 
l2012l l. However, there is growing evidence suggesting that 
the time scale for disc dissipation could depend on stel¬ 
lar mass, because the fraction of objects that retain cir- 
cumstellar discs increases f or those with lower masses (e.g. 
iLuhman fc Mamaiekllioi^ l. Such a dependency could have 
remarkable implications on the efficiency for the formation 
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of giant an d terrestrial planets aroun d stars of different 
masses (e.g. IPascucci et al. I l2009l.l2013ll . 

The first evidences of such a dependency came from 
studies in the stellar mass domain. iHernandez et all (l2005l l 
found that the inner disc frequency at ages between 3 
and 10 Myr in intermediate-mass stars is lower than in 
low-mass stars and suggested that it could be a conse¬ 
quence of a more efficient mechanism of primordial disc 
dispersal in the intermediate-mas s stars. These result s 
were subsequently s uppor ted by Meeeath et al.l LoOSh . 
ISicilia- Aguilar et al. _(|200dl and iLada et al.l ll200^ . Fur¬ 
thermore, Carpenter et al.l (l2006l l studied the number frac¬ 
tion of discs surrounding stars in the wide mass range 
0.1 M/Mq 20 in the Upper Sco association, finding that 
the fraction of primordial optic ally thick discs d e creases 
as the mass of the star increases. IHernandez et all (l2007al l 
studied stars with spectral types K6 to M4 in the Orion OBI 
association and found the maximu m of the disc freq u ency in 
MO. In the younger a Ori cluster IHernandez et al.l (l2007bll 
reported disc fractions of 10 % for M > 2 Mq and 35 % 
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for 0.1 < M/Mq < 1. This tendency is a lso supported by 
result s in the A Orionis association from iHernandez et al] 
(l2010l l who found a disc fraction of ~ 6 % for K-type stars 
and ~ 27 % for M5 stars or later. 


The u nprecedented capabilities of the Spit zer Space 
Telescope jFazio et al.l 200^: Rieke et al.l 120041 ') and the 
WISE survey I Wright et al.ll201ol 'l allowS t o extend these 
studie s down to the sub-stellar mass domainjLuIiman^_et_al 


200^ complemented the results from IHernandez et al 
(l2007lJ ) for the cr Ori cluster finding that 60 % of the BDs 


shows IR excesses consistent with discs. iRiaz fc GizisI ll2008ll 
studied the TW Hya association and found that ~ 60 % of 
the BDs show IR excesses , again st ~ 24 % of the VLMS 
and iRiaz. Lodieu fc GizisI ll2009ll suggest that the longer 
disc lifetimes in TW Hya could be a co nsequence of it s lower 
BD spatial density. In the same region. iMorrow et al.lll2008ll 
showed that the BDs having irradiated accretion discs do 
not show the silicate emission found in d iscs around VLMS 
(lUchida et al.l l2004l : iFurlan et al.l l2007ll and in younger 
VLMS and BDs. They interpreted these results as an indica¬ 
tion that grain growth occurs more rapidly in discs around 
BD than in those around stars or th at grains grow faster 
at sm aller disc radii as suggested by iKessler-Silacci et al.l 
(l2007l ~) and ISicilia-Aguilar et al.l (l2007l ~). 


Parti cularly, the ^ 5 My r old Upper Sco associa¬ 
tion (e.g IPreibisch et al.l I2OO2I ') has been the subject of 
several studi e s abo ut the fraction of BD harboring discs. 
IScholz et al.l (l2007l l studied VLMS and BD with spec¬ 
tral types from M5 to M9, finding a disc frequency of 
~ 37 %. They also reported that ~ 30 % of such objects 
also show Hq em issions consistent with active accretion. 
iRiaz et ^ 1 I 2 OI 2 II compiled all the spectroscopically con¬ 
firmed BDs and identified discs based on WISE photom¬ 
etry. They found a disc frequency of ~ 28 % and that 
half of the VLMS and BD harboring discs also show sig¬ 
natures of accretion. They did not find any dependence 
of the disc fraction with the stellar number density or 
the BD/star number ratio and suggested that the differ¬ 
ences in disc life times could also be a consequence of 
different BD formation mechanisms and/or different ini¬ 
tial disc fractions. Cons idering the latest age estimate fo r 
that region 11 Myr: IPecaut. Mamaiek fc Bubarll20l'^ . 
iLuhman fc MamaiekI (l2012l ~) found that the disc fraction of 


objects with masses 0.01 M/Mq ^ 0.2 reaches ~ 25 %, 
which indicates that such primordial discs co uld survive for 
at lea st ^ 10 Myr. Recently, in this region, IPawson et al.l 
(I 2 OI 3 I I reported that 23 % of the BDs are class H and 
that 19 % are transitional discs. They also compared 
their fractions with those for K and M-type stars in Up¬ 
per Sco, Chal, IC348 and a Ori ( [imhrtian et al.l 20051 : 


Dam janov et al.l l2007l : iLada et al.l I 2 OO 6 I : Hernandez et al.l 
2007al . respectively') and argued that the correlation be¬ 


tween the disc fraction and the stellar mass is not statis¬ 
tically significant, and that the average lifetimes of discs 
around such s tars could not depend on the stellar mass. 


I 2 OIC 1 I: iPetr-Gotzens et al.1 

20111. IRAG-Snitzer and WISE 

nhotometrv Downes et al. 

( 2 OI 4 II studied the 7 Mvr old 


25 Orionis group and Orion OBla finding a number fraction 
of candidates low-mass stars (LMS) showing IR excesses be¬ 


tween ~ 7 % and ~ 10 % while for BDs candidates the 
number fraction increases up to ~ 20 % to ~ 50 %. 

In summary, there are several indications on a possible 
dependence of the characteristic time scale of disc dissipa¬ 
tion with stellar mass, that extends down to the sub-stellar 
mass regime, although some issues related to the estimation 
of ages and the uncertainties in the number fractions need 
to be clarified. These issues could be addressed by study¬ 
ing different disc indicators in slightly more evolved regions 
(t ^ 5 Myr) with consistent age estimations. This, together 
with the use of numerous samples of LMS and BDs and a 
thorough statistical treatment, can provide a robust analy¬ 
sis of the variations of disk fractions with age. 

In this paper we present an optical spectroscopic 
and optical/near-IR photometric study of 15 BDs with 
masses 0.01 ^ M/Mq ^ 0.07 (spectral types between M6 
and M9), together with the sample of 77 LMS with masses 
0.1 M/M(7] ^ 0.8 ( spectr al types between MO.5 to M5.5) 

from I Downes et ahl ll2014ll . These 15 BDs were spectroscop¬ 
ically confirmed as members of the 25 Orionis group and 
its surroundings in Orion OBla, from an initial sample of 
21 BD candidates. These regions have an age of ^ 7 Myr , 
estimated from LMS and BD samples llBriceno et al.ll2005l . 
I 2 OO 7 I : iDownes et al.11^0141) . all of which are consistent. We 
obtain the disc number fractions of BD and LMS from 
two different indicators: IR-excesses and spectroscopic sig¬ 
natures of ongoing accretion. We compute the fractions on 
both sides of the sub-stellar mass limit following the same 
procedure, and provide a statistically robust treatment that 
allows us to compute the probability that the disc fractions 
of LMS and BD are different. 

The paper is organized as follows: In Section [2] we de¬ 
fine the sample and describe the photometric database, the 
spectroscopic observations and data reduction. The mem¬ 
bership diagnoses are discussed in Section [3] In Section 2] 
we classify the new members as class H, evolved discs and 
class HI according to the IR photometric signatures. In Sec¬ 
tion |S] we analyse the spectroscopic signatures of accretion 
and classify the new BDs as Classic T Tauri star (CTTS) or 
Weak T Tauri star (WTTS) sub-stellar analogous. In Sec- 
tion[7]we comment on particular objects and the discussion 
and conclusions are summarized in Section [8] 


2 THE SAMPLE, OBSERVATIONS AND DATA 
REDUCTION 

Eor the present work we studied 21 targets from the sam¬ 
ple of photometric BD candidates with expected spectral 
types between M6 and LI, obtained during the survey 
of the 25 Orionis group and its surroundings carried out 
bv iDownes et al] (|2014| L The candidate selection was per¬ 
formed based on their position in colour-magnitude dia¬ 
grams that combine I-band optical p hotometry from the 
CIDA Deep Survey of Orion fCDSO. iDownes et al.|[2014l l 
and near IR photometry in the J, H and Ks bands from 


the Visible and Infrared Survey Telescope for Astronomy 

(VISTA. lEmerson et al.1 

2004 lEmerson & Sutherlandll2010l: 

iPetr-Gotzens et al.1 2011 

)• 


We have carried out the spectroscopic obs ervations of 
the 2 1 candidates with the OSIRIS instrument llCepa et al.l 
I 2 OOOII at the Gran Telescopic de Canarias (GTC). The 21 
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N [4’x4’] 


0.0 0.2 0.4 0.6 0.8 1.0 1.2 



84.1 83.6 83.1 82.6 82.1 81.6 81.1 80.6 


an 


Figure 1. Spatial distribution of the sample of the new 15 BDs confirmed in this work (labelled white circles) and those confirmed as 
field dwarfs (labelled gray circles). The labelled big asterisks indicate stars o f the Orion belt and the dashed circle roughly encloses the 
~ 5 Myr old Orion OBlb sub-association located at 440 pc from the Sun ijBriceno et al.|[2Q05lb The background gray scale indicates 
the number o f sour ces per 4' x A' bin, of LMS and BD candidates from Downes et al. (in preparation) in Orion OBlb, and from 
IPownes et alJ l|2014^ in the 25 Orionis group (over- density in the uppe r right corner) and its surroundings. The labelled black squares 
indica te the known stellar groups in the area from [Kharchenko et al] ll2005l i together with the 25 Orionis group from iBriceno et al.l 
ll2nn5h . The arrow indicates the displacement that a hypothetical star would have if moving from Orion OBlb toward the 25 Orionis 
group during 5 Myr with a velocity of 2 km s”^ (see Section [3] for discussion). Empty squares indicate the BDs classified as class II or 
evolved discs and the empty circles the CTTS as discussed in Sections [4] and [6] 


candidates analysed in this work were selected as follows: 
First we chose 16 candidates with photometric colours con¬ 
sistent with members of Orion with spectral types M7 to 
LI. These were randomly selected from the full catalog of 
BD candidates, i.e. without imposing any further cuts that 
might bias the sample towards objects harbouring discs. 
Out of these 16 candidates, 5 objects had another M6 to 
M7 candidate at a distance smaller than 7.4', which is close 
enough to be observed simultaneously with the single long- 
slit configuration of OSIRIS. In this way, we obtained the 
spectra for 21 candidates with estimated spectral types be¬ 
tween M6 and LI, in only 16 observations. The spatial dis¬ 
tribution of the sample with respect to the Orion OBI as¬ 
sociation is shown in Figure [T] 

The observations were performed in service mode dur¬ 
ing March, October and December 2012 and October and 
November 2013 as part of the guaranteed Mexican time with 
GTC. We used the OSIRIS spectrograph in the single long- 
slit configuration with a 1" wide slit and the R500R grism, 
which results in a 4800 ^ A/A )$ 10000 wavelength coverage 


with a dispersion of 4.88 Kjpixel and a nominal resolution 
of 587 at 7319 A. Sky flats, dome flats, bias frames and 
several spectra of comparison lamps for wavelength calibra¬ 
tion were obtained each night. The observing log is shown 
in Table [T] 

The spectra were reduced using standard IRAF rou¬ 
tines consisting of bias subtraction, fiat-fielding, instrumen¬ 
tal response correction, spectrum extraction, removal of at¬ 
mospheric spectral features and wavelength calibration. The 
wavelength calibration was performed with a mean accuracy 
of ~ 0.3 A and cosmic rays were successfully removed during 
the extraction of the spectra. 

We computed the spectral types fo l lowing the semi- 
automated scheme of iHernandez et all (|2004||H and the 
equivalent width of the Ha line was measured by a linear 


^ This procedure was performed us¬ 

ing the code SPTCLASS, available at 
http: //www. astro.Isa. umich.edu/~hernandj/SPTclass/sptclass.html 
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Table 1. Observing log for the observations with OSIRIS at GTC 


Target^ 
ID “ 

Beginning of the observation 
[UT] 

^int 

[b] 

Air mass 

Seeing 

["] 

Observing conditions 

1 

2012-03-12X21:02:06.050 

2400 

1.26 

1.5 

clear/dark 

2 

2012-10-08X03:28:22.655 

2600 

1.28 

1.1 

clear/gray 

3 

2012-10-08X04:32:45.348 

2600 

1.15 

1.1 

clear/gray 

4 

2012-12-09X03:15:59.351 

2600 

1.26 

1.0 

spectroscopic / dark 

5 

2012-12-17X01:06:16.672 

2500 

1.12 

1.0 

photometric / dark 

6 

2012-12-17X03:24:48.060 

2500 

1.41 

1.0 

photometric / dark 

7 

2012-12-21X21:25:13.901 

2600 

1.61 

0.8 

spectroscopic / gray 

8 

2013-10-30X06:02:09.175 

2100 

1.29 

0.8 

clear/dark 

9 

2013-10-13X05:02:08.778 

2500 

1.13 

0.7 

clear/dark 

10 

2013-10-13X05:02:08.778 

2500 

1.13 

0.7 

clear/dark 

11 

2013-10-14X05:01:06.925 

2300 

1.12 

0.9 

photometric/bright 

12 

2013-10-14X05:01:06.925 

2300 

1.12 

0.9 

photometric/bright 

13 

2013-11-06X05:05:57.362 

2300 

1.20 

0.7 

spectroscopic / dark 

14 

2013-11-05X04:39:54.909 

2400 

1.15 

0.7 

spectroscopic / dark 

15 

2013-11-05X04:39:54.909 

2400 

1.15 

0.7 

spectroscopic / dark 

16 

2013-11-06X01:50:38.996 

2300 

1.26 

0.8 

spectroscopic / dark 

17 

2013-11-06X01:50:38.996 

2300 

1.26 

0.8 

spectroscopic / dark 

18 

2013-11-06X04:15:15.002 

2300 

1.14 

0.9 

spectroscopic / dark 

19 

2013-11-06X04:15:15.002 

2300 

1.14 

0.9 

spectroscopic / dark 

20 

2013-11-07X02:23:33.637 

2200 

1.19 

0.6 

spectroscopic / dark 

21 

2013-10-30X06:02:09.175 

2100 

1.29 

0.8 

clear/dark 


^ These IDs allow the identification of the targets in the figures and tables of the article. 
Additional designations from literature are included in the electronic version of Table [3] 


Table 2. Spectroscopic catalog for the candidates observed with Osiris. 


Xarget 

ID 

SX 

CaH^ 

VQi 

VO^ 

KI 

Nal 

WHa 

[A] 

Av 

[mag] 

Membership 

1 

M7.5±1.0 

-1 

-1 

-1 

-1 

-1 

0 

0.37 

field 

2 

M9.0±1.0 

1 

1 

1 

0 

1 

-21 

0.02 

member 

3 

M8.0±1.5 

0 

1 

1 

-1 

1 

-381 

0 

member 

4 

M6.0±0.5 

-1 

1 

-1 

1 

1 

-27 

6.29 

member 

5 

MS.Oitl.O 

1 

1 

1 

0 

1 

-27 

0 

member 

6 

M8.5±1.0 

0 

1 

1 

0 

1 

-226 

0 

member 

7 

M9.0±1.5 

-1 

0 

0 

0 

1 

0 

0 

field 

8 

M7.5±2.0 

0 

1 

1 

0 

1 

-18 

0 

member 

9 

M8.0±0.5 

1 

1 

1 

0 

1 

-81 

1.36 

member 

10 

M7.0±0.5 

1 

1 

1 

1 

1 

-10 

0 

member 

11 

M7.0±0.5 

0 

-1 

0 

1 

0 

-17 

1.61 

field 

12 

M7.5±0.5 

1 

1 

1 

0 

1 

-19 

0 

member 

13 

M9.0±0.5 

0 

0 

1 

0 

1 

0 

0 

field 

14 

M8.0±1.5 

-1 

1 

1 

0 

1 

-50 

0 

member 

15 

M6.0±0.5 

1 

0 

1 

1 

1 

-14 

0.61 

member 

16 

M7.0±0.5 

-1 

-1 

-1 

0 

-1 

0 

0 

field 

17 

M7.0±0.5 

1 

1 

1 

1 

1 

-27 

0 

member 

18 

M8.0±0.5 

1 

1 

1 

1 

1 

-44 

0 

member 

19 

M7.0±0.5 

1 

1 

1 

1 

1 

-18 

0.21 

member 

20 

LI.Oil.0 

0 

-1 

0 

-1 

-1 

0 


field 

21 

M7.5i0.5 

1 

1 

1 

1 

1 

-213 

0.12 

member 


^ The flags indicate if the spectral feature is consistent with a young BD (1), with a field 
dwarf (-1) or not conclusive (0). See Section [3| for details. 
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Figure 2. Spectra and spectral features for the sample of the young BDs analysed in this work. Lower pannels show the indicators of 
youth (VO, Nal, KI and Cal) and ongoing magnetospheric accretion (Ha). Solid lines indicate the spectra for the original resolution 
of ~ 9 A (solid gray line) a nd for ^ 16 A res o lution (solid black line). Dashed and dotted lines indica te respec t ively a field dwarf of 
the same spectral type fromlKirkoatric k et al.l ijlQQgh and a young dwarf of the same spectral type from llmhmanl l|2000lL [Briceno et al.l 
ll2002t) . [Luhman et al.l ll2003h and lLuhmanl ll2004h both at a resolution of 16 A. The vertical solid gray lines indicate the wavelength 
range of the corresponding feature. This plot corresponds to the BD number 6. The plots for the remaining BDs confirmed as members 
are available in the electronic version of the article. 


fit of the continuum performed with the splot task from 
IRAF. Figure [2] shows the spe ctra of the new member s 
compare d with standards from Kirkpatrick et all d 19991) . 
Luhmanl (l2OQ0l ) .iBriceno et al.l l|2002h . lLuhmanl ~( 20041 ) and 


Luhman et alJ (l2005l 'l and Table [2] shows the derived spec¬ 


tral types, visual extinctions, equivalent widths of the Ho? 
emission line, the surface gravity indicators and the final 
membership diagnosis we will explain in Section [S] 

In order to detect possible IR-excesses from the spec¬ 


tral energy distribution (SED), we complemented the I- 
band photometry from CDSO and the J, Z, Y, H and 
Ks-band photometry from VISTA with additional pho¬ 
tometry at 3.6, 4. 5, 5.8 and 8.0 um from IRAC-Spitzer 
observations from iHernandez et all ll2007ah and Briceno 
et al. (in preparation), and at 3.4, 4.6, 12 and 22 /rm 
from the WISE All-Sky Source Catalog l|Wright et al.l 
l20ld l. Additional photometric data in the g, r, i and 
z-bands from the Sloan Digital Sky Survey Catalog 
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Data Release 8 llAdelman-McCarthv fc et alJ I 2 OIIII were 
obtained using; the VizieR virtual observatory system 
(lOchsenbein, Bauer fc MarcoutI [ 2 OOC 1 I I. Because of the dif¬ 
ferent sensitivity limits and spatial coverage of each survey, 
not all the photometric bands are available for all the can¬ 
didates. The available photometric information for the new 
BDs is summarized in Table [3] 


As discussed by Ijarrett et al.l (1201 il l there is an off¬ 
set between WISE 3.6 /rm and IRAC 3.4 /im magnitudes, 
as well as between WISE 4.5 fj,m and IRAC 4.6 fim mag¬ 
nitudes. Both offsets occur approximately for ms.e > 
14 magnitudes and m 4.5 > 13 magnitudes in such 

a way that the magnitudes from WISE are increasingly 
fainter than t he co rresponding from Spitzer. According to 
Ijarrett et al.l 1 I 2 OIIII these biases occur as a consequence of 
an overestimation in the background levels of the WISE im¬ 
ages. We found exactly the same behaviour affecting all of 
our faint objects with available photometry from WISE and 
Spitzer. In order to correct the photometry from WISE, we 
selected 10000 objects in the Spitzer helds with avail¬ 
able photometry from WISE and fit the residuals of the 
corresponding magnitudes as a function of one of them. Af¬ 
ter that, we applied the corresponding offset to the WISE 
photometry as a correction of the reported magnitudes. Be¬ 
cause IRAC-Spitzer has no measurements around 12 /rm 
it was not possible to establish if the background in the 
WISE 12 ^m pass-band is also underestimated. Only two 
BDs (numbers 4 and 18) have measurements in the WISE 
12 /rm pass-band with reasonable SNR. We disregard all 
the measurements in the WISE 22 /rm pass-bands because 
of the very low SNR. 


3 MEMBERSHIP DIAGNOSES 


Young BDs are still contracting and their surface grav¬ 
ity is lower than for an old held dwarf of the same ef¬ 
fective temperature. Since the contaminatio n of the can¬ 
didat e sample comes from old held dwarfs dPownes et al.l 
I 2 OI 4 II . we evaluate the membership of the candidates to 
the 25 Orionis group or Orion OB la, according to sev¬ 
eral s£ectralMeati^s sensitive to surface gravity follow¬ 
ing iMcGovern et ah] (l2004|j . We considered the following 
spectral features: (z) The CaH (A6750 to A7050) molecu¬ 
lar band, which is stronger in held dwarfs; (ii) the VO^ 
and VO^ molecular bands (A7300 to A7500 and A7800 to 
A8000 respectively), which are weaker in held dwarfs; (Hi) 
the atomic lines KI A7665, 7699 and Nal A8183, 8195, which 
are expected to be stronger in held dwarfs. Other spectral 
features sensitive to surface gravity, such as the absorption 
lines Rbl A7800, 7948 and Csl A8521, were marginally de¬ 
tected only in some of the spectra because of the low SNR 
and wavelength resolution and did not allow for a reliable 
estimation of the surface gravity, hence they were not con¬ 
sidered further. 

These 5 spectral features were compared in th e spectra 
of our c andidates to t h ose o f young dwarfs fr om LuhmanI 
(l2000l~). Briceno et al.l (l2002ll . iLuhman et al.l (l2003l~) arid 
iLuhman ( 2004h . and held dwarfs from iKirkoatrick et al.l 
(I 1999 I I of the same spectral types. We found a very good 


general agreement between the different surface gravity in¬ 
dicators as we show in Table [ 2 ] and Figure [21 Table [2] sum¬ 
marizes, for each of the 21 candidates, whether each of the 
5 spectral features is consistent with a young or held BD of 
the corresponding spectral type. We classihed as members 
those candidates showing at least 3 spectral features consis¬ 
tent with low surface gravity. Using these criteria we have 
conhrmed 15 members from the sample of 21 photometric 
candidates. 

We stress that such a selection is also supported by: 
(i) all the objects classihed as young BDs show Ha line in 
emission and {ii) their visual extinction is consistent with 
the mean value for kno wn members in 25 Ori and Orion 
OBla (Ay ~ 0.4 ± 0.3 iDownes et al. I I 2 OI 4 I : iBricefio et al.l 
I 2 OO 7 I . I 2 OO 5 I I. The visual extinction for each candidate was 
obtained from the observed I-J colour and the intrinsic I- 
J colour that corresponds to the spectral type. We used 
the intrinsic co lo urs to spectral type r elatio nships used by 
iLujhmar] (Il999l l. iBriceho et al.l ll2C)02h and Luh man et aP 
( 20031) . designed to match the iBaraffe et"^ ( 1998h tracks 
such that the components of G G Tau appear c o eval o n the 
H-R diagram, as explained by Luhman et al. (l2003h. We 
used the extinction law from iGardelli. Clayton fc Mathis! 
(Il989f) assuming Rv = 3.09. The resulting extinctions are 
shown in Table [2] and are consistent, within the uncertain¬ 
ties, with previous determinations. The only exceptions are 
for the BDs 4 and 9 and are discussed in Section [T] 

The selected young BDs could be members of the 25 
Orionis group or other regions of Ori on OBla which es¬ 
sentially have the s ame age (~ 7 Mvr. iBriceno et 111 ! I I 2 OO 5 I: 
iDownes et al. Ilioli). Nevertheless, the BDs could also be 
members of a different sub-region within the Orion OBI as¬ 
sociation. Particularly we are interested in determining if 
some of them could belong to the younger and relatively 
close (~ 3°, FigurelD) Orion OBlb (~ 5 Mvr. IBricefio et all 
I 2 OO 5 I) . However, we conclude the contamination due to 
members of the nearby younger Orion OBlb must be very 
low for the following reasons: 

(i) The spatial distribution of the new BDs is con¬ 
sistent with those from previously known members of the 
25 Orionis group and Orion OBla of earlier spectral types 
([Briceno et al.l[2005l: IDownes et al.l[^014l: iKharchenko et ah! 

l2005f ). As shown in the spatial overdensities of Figure [T] the 
candidate LMS and BDs of Orion OBla and the 25 Orionis 
group are spatially distinct from the region populated by 
Orion OBlb sources. 

(ii) Considering a velocity dispersion of ^ 2 km s“^, 
the angular separation between the sub-association Orion 
OBla and OBlb and the i r dist ances to the Sun 360 pc 
for OBla (IBricefio et al.l l2007l) and 440 pc for OBlb 
dBriceno et al.| 20051)), a me mber of the younger OBlb 

5 Mvr. IBricefio et eJ would need more than 5 Myr 

to escape and end up in the line of sight to the 25 Orionis 
group. Thus, the possible contaminants could, at most, af¬ 
fect the eastern part of the survey. However, if the contam¬ 
ination from OBlb were significant an overdensity of class 
II, evolved discs and/or accretors (see Sections |4] and [6|) 
would be expected in the eastern part of the survey closer 
to Orion OBlb, which is not observed in Figure [T] There¬ 
fore, we consider it much more likely that these young BDs 
are members of the 25 Ori group or the Orion OBla sub¬ 
association, than Orion OBlb contaminants. 
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Table 3. Photometric catalog of the new BDs confirmed as members of 25 Ori or Orion OBla 


ID_^ 

RA 

DEC 

I 

J 

Z 

K 

3.6 fim 

4.5 fim 

2 

80.791100 

1.714128 

20.79 

17.49 

19.77 

16.42 

16.27 

16.72 

3 

80.850308 

1.790936 

19.93 

17.26 

19.41 

16.21 

15.68 

15.00 

4 

80.894376 

1.846516 

19.79 

16.52 

18.94 

15.12 

14.80 

14.31 

5 

80.951425 

1.809198 

20.20 

17.30 

19.49 

16.23 

15.98 

15.62 

6 

81.243647 

1.733360 

20.32 

17.35 

19.49 

16.24 

16.07 

15.65 

8 

81.514956 

0.737391 

19.25 

16.67 

18.39 

15.84 

15.43 

15.20 

9 

82.375058 

0.790957 

20.76 

17.42 

19.58 

16.32 

16.02 

15.71 

10 

82.307508 

0.852826 

18.45 

16.13 

17.83 

15.26 

15.09 

14.88 

12 

81.956878 

1.205257 

18.35 

16.21 

17.70 

15.38 

15.14 

15.04 

14 

82.675353 

1.859462 

20.21 

17.41 

19.57 

16.33 

16.07 

15.71 

15 

82.647725 

1.775001 

18.15 

16.02 

17.41 

15.22 

14.98 

14.77 

17 

81.909127 

0.961591 

18.26 

15.99 

17.77 

15.09 

14.72 

14.02 

18 

81.673711 

0.637468 

19.75 

17.04 

19.18 

16.03 

15.51 

15.02 

19 

81.648441 

0.600469 

18.70 

16.22 

17.96 

15.38 

15.10 

14.83 

21 

81.549774 

0.820961 

19.46 

16.67 

18.52 

15.73 

15.51 

14.89 


“ The complete version of the table is available in the electronic version of the article. 


(Hi) A particular case is the new BDs 8, 18, 19 and 21 
which are close to the ste ll ar gro up ASCC18 from the cat¬ 
alog of iKharchenko et al.l (|20n5l ~l who report it as a group 
in the Orion Q Bl association. On t he ba sis of one interme¬ 
diate mass star, IKharchenko et ahl (l2005l l estimated an age 
of 13 Myr for ASCC18 which is slightly older than the 25 
Orionis group and Orion OBla, something which does not 
affect the conclusions presented here. 

The position of the new BDs in the H-R diagram is 
shown in Fig ure El together with the LMS of the 25 Orionis 
group from ll Downes et al. I [2^. We computed the effec¬ 
tive tem^eraturesbyh^rpolation of the spectral types into 
the iLuhman et al.l 1 I 2 OO 3 II relationships for young BDs and 
we computed the bolometric luminosities from the dered- 
dened I-band magn itudes, assuming a distance of 360 pc 
(iBriceno et al.ll2007l l and the bolometric corrections from 
IPahn et al.l l ^02ll . The masses were derived by interpo¬ 
lation of the luminosities and the effective tempera tures 
into the DUSTY models from IChabrier et ahl (I 2 OOOI I and 
are shown in Table [5] 

Finally, based on si mulations perform ed with the Be- 
sangon Galact i c mo del (iRobin et al.ll200^ . we showed in 
iDownes et al.l (l2014l l that in the spectral-type range con¬ 
sidered here, our photometric selection procedure for can¬ 
didates in the 25 Orionis group and its surroundings has 
a mean total contamination of 25 %, composed solely 
of foreground field dwarfs. In this work we spectroscopi¬ 
cally confirmed that 71.4j^g Q % of the candidates are real 
members and the remaining targets show spectral features 
consistent with the high surface gravity expected from old 
dwarf stars from the field, which perfectly matches the ex¬ 
pected contamination fraction. 


4 PHOTOMETRIC SIGNATURES OF DISCS 

It is well established that some youn g BDs show IR excesses 
indicative of circumstelar discs (e.e. iLubmanll^PlJ l . In this 
section we classify the discs surrounding the new BDs ac¬ 


cording to their excesses in the wavelength range A > 2 /rm. 
For the classification of the discs we assume the general 
scheme for their evolution around LMS and VLMS (e.g. 
iLada et al.ll200d ~l in which they evolve from class 0 to I, then 
to class II, then passing through the evolved disc stage (e.g. 
iHernandez et al]l2nn7al L perhaps passing t hrough the pre- 
trans itional to the transitional stages (e.g. lEsoaillat et al.l 
l2007l 'l depending on the mass of the disc, and ending as class 
III objects. Thus, we follow the same classification scheme 
in order to compare the evolution of discs around objects 
at both sides of the sub-stellar mass limit. 

The evolutionary stage of the discs can be inferred from 
the slopes of the spectral energy distributions (SED) at dif¬ 
ferent wavelength ranges or equivalently from its distribu¬ 
tion in selected colour-colour diagrams. We performed the 
classification according to the SEDs of each new member 
and plotted the results in a set of selected colour-colour di¬ 
agrams in order to show the distribution of the BDs and 
how it supports their classification through the SEDs. All 
the SEDs include measurements at the I, Z, Y, J, H, Ks, 
3.4 /rni and 4.6 fj,m photometric-bands, although the differ¬ 
ences in sensitivity and spatial coverage of the WISE, IRAC 
and SDSS surveys give us information only for a subset of 
the complete sample in the remaining photometric-bands 
mentioned in Section [2] 

We consider as IR excesses all the SED points 
for A ^ 2 /rm showing fluxes above the photospheric 
emissions predicted by the BT-Dusty model from 
lAllard. Homeier fc Frevtad (|2012h which was developed for 
near-IR studies of BDs with Te// > 1700A'. Our procedure 
was as follows: First, we estimate the photospheric emis¬ 
sion for the BDs by fitting their extinction corrected SEDs 
to the BT-Dusty model. For all the fits we used the I, Z, 
Y, J and H band-passes which were complemented with 
the g, r, i and z photometry when available. In this way, 
we perform the fits with a minimum of 5 and a maximum 
of 10 band-passes in a wavelength range where the IR ex¬ 
cesses are not expected to occur and where the peaks of the 
photospheric contribution to the SEDs are expected. Before 
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Figure 3. H-R diagram of the new young BDs confirmed in this work (labelled circles) together with LMS from lPownes et alJ ll2014h 
(asterisks). The big circles indicate the BDs placed in the 25 Orionis group. The evolutionary tracks for 0.8, 0.7, 0.6, 0.5, 0.4, 0.3, 0.2, 
0.1, 0.072 (dashed), 0 .06, 0.05, 0.04, 0.03, 0.02 and 0.01 Mq, and isochrones for 1, 3, 5, 10 and 30 Myr from the DUSTY models of 
IChabrier et alj ll2000l~) are also indicated. We applied a fiducial 40 K offset in temperature to the BDs 9 and 14 for clarity. 


the fits, the SEDs were corrected by reddening using the 
Av obta ined as explained in Section [3] and the extinction 
law fr om iFitzpatrickI (1 19991 ) improved bv lindebetouw et al.l 
(|200 !tI ) in the infrared (A > 1.2 pm), with the parameter 
R=3.02. We make use of the Virt ual Observatory SE P An¬ 
alyzer (VOSA) fitting tools from iBavo et al.l (l2008l) in or¬ 
der to apply the reddening corrections and make the fits 
to the lAllard. Homeier fc Frevtad (l2012l ) models. The effec¬ 
tive temperatures obtained from the fits are consistent with 
tho se obtained by t he int erpolation of the spectral types in 
the lLuhman et all l|2003h relationships. 

We classify as class II those BDs showing excesses at the 
K-band or longer wa velengths, consisten t with the median 
SED for Taurus from iFurlan et all ll2006t) . We considered as 
evolved discs those BDs showing excesses for wavelengths 
longer than the Ks-band but clearly weaker than the ex¬ 
cesses observed for class II. Finally we classify as BDs of 
class III those showing a purely photospheric SED consis - 
tent with the fits to the lAllard. Homeier fc Frevtad l|2012l ) 
models. 

With the available photometry we cannot reliably de- 
tect BDs surrounded b y transitional or pre-transitional discs 
llEspaillat et al.l [20071 ) because the spectral coverage in the 
wavelength range 8 ^ A/pm .$11 is not homogeneous for 


most of the sources. Figured shows the SEDs for all the new 
confirmed BDs, toge ther with the resulting fit to the p hoto- 
spheric models from I Allard. Homeierj^^jTey tag j[2^2l) and 
the median SED for Taurus from iFurlan et al.l ( 20061 ). 


The IR excesses are also detectable in colour-colour 
diagrams supporting the results obtained from the SEDs. 
Figures [5] and [6] show a selection of colour -colour diagrams 
some of those showing the discs locii from Hartmann et ^ 


(I 2 OO 5 I ) , iLuhman et al.l (I200R| ) and iLuhman et al.l 1 201(t l. 


Summarizing the results from our classification, out of 
the 15 new BDs, we have found 3 objects showing IR ex¬ 
cesses consistent with BDs of class II, 8 consistent with 
evolved discs and 4 showing purely photospheric SEDs 
that were cla s sified as BDs of class III. Using data from 
iDownes et al. IHli), we have followed the same procedure 
to classify IR excesses observed in a sample of 77 LMS 
with spectral types between MO.5 and M5.5 and masses 
0.1 $ M/Mq $ 0.8, confirmed as members of the 25 Orionis 
group. We find 3 of class II, 11 having evolved discs and 63 
of class III. These results are summarized in Table [d] 


An important issue is that our selection of the photo¬ 
metric candidates was performed inside a set of locii in the 
colour-magnitude diagrams I vs. I-J, I vs. I-H and I vs. I-K 
that enclose all the LMS and BD expected to have purely 
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Figure 4. The black dots and the black solid lines indicate the observed SEDs of t he new BDs. The gray lines show the photospheric 
SED that results from the interpolation of the fluxes in the BT-D usty models from lAllard. Homeier fc FrevtaS ll2012h and the dotted 
lines indicate the median SED for Taurus from iFurlan et a i] ||2006f) . All the SEDs are normalized to the flux at the J-band. Each pannel 
shows the ID of the BD and its classification as class II, evolved disc or class III. The BD number 4 shows clear excesses for H-band 

and longer wavelengths and th e nu mber of points at shorter wavelengths does not allow for a reliable fit into the photospheric models. 
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Table 4. Number fraction of BD and LMS classified as class II, evolved disc, class III, WTTS and CTTS. 


Classification 

BD ( 

<0.072 Mq) 

LMS (0.1 

<M/M0<O.8) 

P>0.1 


[N] 

[%] 

[N] 

[%] 


Class II 

3 

2o.oj:®® 

3 

3-9t?:6 

0.9908 

Evolved 

8 

CO q-|-10.9 
oo.o_ii 0 

11 

14.313-8 

0.9987 

CTTS 

5 

qq q-|-10.8 

00.0_g g 

3 

0 9+2.4 

0.9993 

CTTS & Class II or Evolved 

5 

qq q-|-10.8 

00.0_g g 

3 

3-9i?:6 

0.9993 

Class III 

4 

26.7^1^ 

63 

81.813® 

0.9988 

WTTS 

10 


74 


0.9809 


S ihotospheric emissio ns as well as those showing IR excesses 
Downes et al.l[20l3l . Then, our procedure do not bias the 
selection of candidates towards LMS or BDs harbouring 
discs, even for those candidates that are fainter than th e 
completeness limits in such diagrams llDownes et alJlioidh . 


= ( 20 . 0 ji 7 g) %, which correspond to the most probable 
values and the respective la confidence intervals, computed 
as the 16th and 84th percentiles of the corresponding PDFs, 
respectively shown with dashed and dotted lines in the up¬ 
per left panel of Figured 


5 COMPUTATION OF DISC FRACTIONS 

In order to compute discs fractions and their uncertainties 
given the limited number of objects with discs in our sam¬ 
ples, we use Bayesian statistics. This allows us to state the 
problem in a general way providing us with a probability 
function for the disc fractions and, more importantly, with 
a quantitative way of estimating the probability that disc 
fractions from two independent samples differ. 

We can think of the disc fraction fd as the probabil¬ 
ity for any single star to harbour a disc. Given a sample 
with a disc fraction fd and a total number of Nt stars, the 
likelihood P{Nd, NT\fd) of observing Nd stars with discs 
is then simply given by the Binomial distribution. Assum¬ 
ing a uniform prior probability distribution for the disc 
fractio n in the range 0 ^ f d ^ 1, from the Bayes’s the¬ 
orem dSivia fc Skillind l200d ') we have the Posterior PDF 
expressed simply as: 


The next question is what is the probability that these 
two populations have differing class II fractions. It is clear 
that in this case the obtained disc fractions are quite dif¬ 
ferent, even taking into account the reported uncertainties, 
but the posterior PDF of Eq. [5] allows us to quantify this 
probability in a general manner, which will be useful for 
less clear cases. Marginalizing the 2D posterior PDF, we 
compute the probability that the class II fractions of 
the two populations differ by more than a certain fractional 
margin x as: 


P>. = 1 - / dXPifJjS) = il+X)fX\{Nd,NT}^’^){3) 


In this case, choosing a fractional difference of 0.1, we find 
P>o.i = 0.9908, i.e. there is a 99.08 % probability that the 
class II fractions of LMS and BD are different by more than 


P{D\Nd,NT)=CfX^{l-fd) 




( 1 ) 


where C is a normalization constant such that 

/ dfdP{fd\Nd,NT) = 1. 

We can now generalize this to the case where we have 
two independent samples, LMS and BDs, and express the 
full posterior probability of fj and ff as the product of 


the two independent probabilities: 

P{fXDi\NX,NX,NX,NX) = P{fX\NX ,NX)P{fX\NX,NX){2 

where NJ and A® are the observed numbers of LMS and 
BDs harbouring discs, in samples with a total number of NX 
and Ay of LMS and BDs respectively. Each of the terms in 
the right hand side in this equation are thus given by Eq. 

m 


10 %. 

Similarly, we use Eq. [5] to compute the Posterior PDF 
for the fractions of evolved discs and class III objects in 
the LMS and BD mass regimes. For evolved discs we find 
fractions of Id.Slg ® % and SS.Slj’j’ Q % for LMS and BDs 
respectively, with a probability P>o.i = 0.9987 (from Eq. 
H that these two fractions differ by more than 10 % (Fig¬ 
ure 0 . For class III we obtain fractions of 8 I. 8 I 4 Q % and 
for LMS and BDs respectively, with a probabil- 
ty P>o.i = 0.9988. The number fractions and Pq.i probabil¬ 
ities are summarized in Table 0 We conclude that the num¬ 
ber fractions of class II, evolved discs and class III clearly 
differ at both sides of the sub-stellar mass limit. This dif¬ 
ferences are statistically robust as supported by the P>o.i 
probabilities. 


Applying this for the LMS and BD samples detailed in 
the previous section, for class II objects we have NJ = 3, 
A® = 3, Ay = 77 and Ay =15 which results in the Pos¬ 
terior PDFs shown in the upper panels of Figure 0 The 
right panel shows isocontours of the 2D posterior of Eq. 0 
The left panel shows the marginal ID Posterior PDFs for 
the class II fraction of LMS (gray solid line) and BDs (black 
solid line). The resulting fractions are fX = ( 3 . 9 l 4 'g) % and 


Finally, we e mphasize that we con firmed the number 
fractions found by iDownes et al.l ll2014h in the 25 Orionis 
group based only on photometric BD candidates, with a 
sample of spectroscopically confirmed members. The frac¬ 
tions we obtained following the Bayesian proc edure for LMS 
are al s o consistent with th ose reported by Bri ceno et ^ 
(l2005l ~). [Briceno et al.l (l2007l i and iDownes et alU 20141 '). The 
final disc classification of each BD is shown in Table [S] 
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Figure 5. Distribution o f the new BDs (labelled circles) and 
LMS from iDownes et alJ ll2014l') (asterisks) in colour-colour di¬ 
agrams. The upper panel shows the observed J-H vs. H - K dia ¬ 
gram includin g thre e phot ospheric locii: f rom [ Luhmanl i|l999l~) . 
iBriceno et alj <|200^ ) and iLuhman et al.l <l20o3^ (dotted line); 
from Ibubman e^^ [ ll201(t ~ (dashed line with labeled squares 
indicating spectral types) and from IPecaut &: MamaiekI ll2013l ) 
(solid line). The dash-dotted li ne ind icates the CTTS locus from 
iMever, Calvet fc Hillenbrandl lll997l for MO stars. In the lower 
panel the shadowed polygon and the dashed polygon indicate 
respectively disc and trans itional disc loc i comp uted from the 
colour intervals defined by iLuhman et alJ ll2niflh . Only BDs 3 
and 4 have information in the 8 pm passband. 


6 SPECTROSCOPIC SIGNATURES OF 
ACCRETION 

The evidence of discs around BDs comes not only 
from IR excesses but also from signatures related 
to magnetospheric accretion such as broad permitted 


Figure 6. Distribution of BDs in colour-colour diagrams. Sym¬ 
bols are as in Figure [5l The shad owed areas indicate the CTTS 
locus from iHa rtmann et al.|l|2005 ) (up per panel), and the excess 
region defined bv iLuhiuarTetaLr i 200fill (lower panel). Only BDs 
2, 3 and 4 have information for the 5.8 pm passband. 


emis sion line profiles and optical continuum veiling 

(e.g. [w hite fc Basrill200,3l:l,Tavawardhan a. Mohantv fc Basril 


I 200 T Muzerolle et al.l 120051 : Luhmanl 20121 ). The number 

fraction of BDs classified as CTTS sub-stellar analogous 
is also an indication of the evolutionary stage of its popu¬ 
lation in the sense that this number fraction is expected to 
decrease with time. In this section we study the equivalent 
width of Ha as a function of the spectral type in order to 
classify the new BDs as sub-stellar analogous of the WTTS 
or CTTS. Also, we compare these results with those from 
the IR excesses indicative of discs discussed in Section |4] as 
well as with stellar counterparts. 

We classified the new BDs as sub-stellar analogous 
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Figure 7. Probability distributions for the fractions of members 
of class II (upper pannels) and evolved disc (lower pannels). The 
left pannels indicate the probability distributions for the frac¬ 
tions of BDs (solid black curve) and LMS (solid gray curve). The 
dotted vertical lines indicate the Icr confidence intervals. The 
right pannels indicate the probability distribution (gray scale) 
as a function of the fractions of BD and LMS. See text for the 
interpretation of the probabilities. 


of the WTTS or CTTS based on the equivalent width 
of the Ha line in emission and their spectral types, 
according to the empirical c l assific ation scheme from 
iBarrado V Navascues Sz Martini (|2003l) in which a BD is a 
CTTS analogue if its Ha emission is stronger than is ex¬ 
pected from purely chromospheric activity which depends 
on their spectral type (see Figure 0 ). There are several op¬ 
tical spectral features which indicate ongoing accretion in 
LMS and BDs, such as H/?, H 7 , Hel A5876 lines in emis¬ 
sion as well as the veiling of the photospheric absorption 
lines dMuzerolle et alj|2003l . [ioOSl i. However, the low SNR 
in the blue-most range of our spectra does not allow for 
a reliable detection of most of these additional lines or to 
measure veiling. The equivalent widths of Ha are shown 
in Table [2] and the corresponding close up of the Ha line 
for each BD in Figure 0 Figure 0 shows the Ha equivalent 
width as a function of the spectral type for the new BDs 
as well as the limit between WTTS an d CTTS proposed by 
IBarrado v Navascues fc Martini (l2003ll . We find that 5 of 
the new BDs show signatures of active accretion consistent 
with CTTS and that the remaining 10 BDs of the sample 
show low Ha emissions as expected in WTTS. Following 
the procedure explained in Section 0 these corresponds re¬ 
spectively to CTTS and WTTS fractions of % and 

66 . 7 jl®g ^9 % among BDs (see Table jH) . _ 

We classified the 77 LMS from (Pownes et ^ lii^ as 
CTTS or WTTS following the same procedure we used for 
BDs. We find 3 CTTS (spectral types Ml, Ml.5 and M5.5) 
and 74 WTTS (spectral types between MO.5 and M5.5) cor- 
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Figure 8. Equivalent width of Ha versus spectral type for the 
BDs confirmed as members. The shadowed area indicates the 
Ho emis sions above the maximum expected fr om chromospheric 
activity llBarrado v Navascues fc MartirJl200.'jh . 


Table 5. Summary of properties for the new confirmed BDs. 


ID 

Log(L/L0) 

Te// 

[K] 

Mass“ 

[MoT 

TTS 

Classification 

2 

-2.7808 

2400 

0.01 

WTTS 

Evolved 

3 

-2.5400 

2710 

0.03 

CTTS 

Class H 

4 

-1.7326 

2990 

0.08 

CTTS 

Class H 

5 

-2.6480 

2710 

0.03 

WTTS 

Evolved 

6 

-2.6024 

2577 

0.02 

CTTS 

Evolved 

8 

-2.4048 

2808 

0.04 

WTTS 

Evolved 

9 

-2.6476 

2710 

0.03 

CTTS 

Evolved 

10 

-2.2248 

2880 

0.06 

WTTS 

Class HI 

12 

-2.0460 

2808 

0.03 

WTTS 

Class HI 

14 

-2.6484 

2710 

0.03 

WTTS 

Evolved 

15 

-2.1845 

2990 

0.08 

WTTS 

Class HI 

17 

-2.1508 

2880 

0.05 

WTTS 

Evolved 

18 

-2.4640 

2710 

0.03 

WTTS 

Class H 

19 

-2.3276 

2880 

0.06 

WTTS 

Class HI 

21 

-2.4896 

2808 

0.04 

CTTS 

Evolved 


“ Masses were estimated from the interpolation of the _ 

Tgff and luminosities into the DUSTY models from IChabrier et al.l 1 I 2 OOOII . 


responding respectively to 3.9jl^ g % and 96.11,'^ % which 
are consistent with the r esults from iBriceno et al.l (120051 1 
and I Downes et al.l ll2014h in the same region. As we did 
for the IR excesses, in Figure 0 we show the probability 
distribution of the fraction of CTTS in the LMS and BD 
mass regimes. We found that within a 10 % margin there 
are, respectively, probabilities of 0.9993 and 0.9809 that the 
fractions of CTTS and WTTS for BDs have different val¬ 
ues than for LMS. Again, these probabilities show that the 
difference observed in the CTTS and WTTS fractions ob- 
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Figure 9. Same as Figure [7] for the fractions of BD and LMS 
classified as CTTS. 

tained for LMS and BDs has a high statistical significance. 

We emphasize that from the 5 BDs showing accretion 
signatures, 2 were classified as class II and 3 as evolved discs. 
In the LMS regime from the 3 objects classified as CTTS, 
2 show IR excesses indicative of class II and 1 was classi¬ 
fied as an evolved disc. Even if we considered as accretors 
only those BDs that were classified as CTTS and class II, 
the fractions at both sides of the sub-stellar mass limit re¬ 
main different: ® % for BDs, 2.6t^3 ^ fo'" LMSs and 

with a probability Pq.i = 0.9833 that both fractions are 
different within 10 %. Finally, the BD 4 was classified as 
CTTS although it is close to th e WTTS/CTTS limit from 
iBarrado v Navascues fc MartfnI (l2003l 'l. Even if we consid¬ 
ered as CTTS only the BDs showing strong Ha in emission 
(numbers 3, 6, 9 and 21), the CTTS fraction is 26.7tlff 
which is still higher than the fraction obtained in the LMS 
regime. 


7 COMMENTS ON PARTICULAR OBJECTS 

(i) The BDs 4, 10, 12, 15, 17 and 19 have information in the 
J, H and Ks-bands from the surveys VISTA and 2MASS. 
From these, the BDs 19 and 17 show evidence of Icr vari¬ 
ability in the three band-passes but such level of variability 
does not chang e the general results present e d her e because 
the fits to the lAllard. Homeier fc Frevtad (l2012fl models, 
changing the photometry randomly within the variations 
we found, do not change significantly. 

(ii) The BD 4 shows an unexpectedly high extinction 
Av = 6.2, considering that the mean extinction towards 
the 25 Orionis group and Orion OBla is Ay ~ 0.3. This 
BD shows a SED clearly consistent with the median SED for 
Taurus and was classified as a class II although the available 
photometry at wavelengths shorter than H-band, where the 
excesses are not detected, is not enough for a reliable fit 
into the photospheric models. Additionally it shows Ha line 
in emission which is consistent with ongoing accretion. We 
speculate that its strong extinction could be a result of the 
disk being edge-on. 

(Hi) The BD 9 shows an extinction Ay ~ 1.36 which 
is slightly higher but it is consistent with mean extinction 
towards the region if the uncertainty in the spectral type is 
considered. 

(iv) The BD 4 was classified as CTTS and the BDs 14 


and 18 as WTTS but there are close to the CTTS/WTTS 
limit proposed by IBarrado v Navascues fc MartfnI (l2003ll . 
New spectroscopic observations with a higher SNR are 
needed to improve their spectral classification (particularly 
BD 14) as well as the measurement of the Ha equivalent 
widths in order to improve their classification as CTTS or 
WTTS sub-stellar analogous. 


8 SUMMARY AND CONCLUSIONS 

We have studied 21 candidate BDs belonging to the 25 Orio¬ 
nis group and the Orion OBla subassociation and spectro¬ 
scopically confirmed 15 of them as new sub-stellar mem¬ 
bers with spectral types between M6 and M9. Compar¬ 
ing the SEDs of the new members at wavelengths beyond 
2 fj,m with those from the BT-D usty photospheric models 
(lAllard. Homeier fc Frevtaall20ld f for the temperatures cor¬ 
responding to their spectral types and with the mean SED 
for Taurus representative of class II (iFurlan et al.l [200^ . 
we detect IR excesses indicative of class II, evolved disc 
and class III. We also classified the new BDs as CTTS and 
WTTS sub- stellar analogous according to th e empi ric limit 
proposed by IBarrado v Navascues fc MartfnI ll2003ll . 

We have applyed a Bayesian analysis to compute the 
number fractions of BDs showing disc and magnetospheric 
accretion signatures and compare them with those for a 
sample of 77 LMS b elonging to the 25 Orionis group from 
iDownes et al.l ll2014l L classified following the same proce¬ 
dures. We find a significant higher fraction of CTTS, evolved 
discs and class II objects in the sub-stellar mass regime 
than in the LMS regime and a fraction of WTTS and class 
III which is lower for BDs than for LMS. All the differ¬ 
ences were found to have a high statistical significance. Our 
main result is that the number fractions of BDs classified as 
CTTS and Class II or evolved disc result to be 33.3tg*g® % 
for BDs while for the LMS is 3.9l^'g %. Such difference has 
a probability P<o.i = 0.9993 of being real. 

The differences between the fractions of discs and ac¬ 
cretion signatures at both sides of the sub-stellar mass limit 
could be interpreted, at least, considering the following 
three scenarios: 

(i) If the disc fraction were not dependent on the mass 
of the central object, a way to mimic the observed fractions 
would be for the formation of BDs to occur over a longer pe¬ 
riod of time than for the LMS. In this scenario, the younger 
BDs would still retain the discs, increasing the mean disc 
fractions for BDs with respect to the LMS. Nevertheless, we 
found ~ 33 % of BDs are CTTS and class II or evolved disc 
which is similar to the fraction observed for LMS at ages 
of ~ 3 Myr (e.g. iHernandez et al1l2007bl . cr Ori; ~ 35 %). 
Therefore, the observed population should have started to 
form ~ 7 Myr ago and extended only for BDs until ~ 4 Myr 
ago. Such an age difference should be reflected in the H-R 
diagram as a larger scatter of the BD locus relative to the 
LMS, which is not observed in the results presented here. 
In addition, extinction patterns commonly found in regions 
of about ~ 3 Myr, such as a Ori or Orion OBlb , are not 
observed in the surveyed area llDownes et al.|[2014l L 

(ii) If the time scale for disc dissipation were the same 
for BDs and LMS, the initial fraction of discs wo uld have to 
be higher for BDs than for LMS, as suggested bv iRiaz et al.l 
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(l2012ll . in order to reproduce the larger disc fractions ob¬ 
served for BDs relative to LMS at a given age. Since we 
are measuring LMS and BD disc fractions at a single age 
of ~ 7 Myr, our observations alone cannot rule out a differ¬ 
ence in the initial disc fractions. However , comparing with 
observa tions in the ~ 3 Myr ol d a Ori fromiH ernande z et al.l 
ll2007bl . ~ 35 % for LMS) and lLuhman et ahT i 20081 . ~ 60 % 
for BDs), we see that for LMS the fraction drops from 
~ 35 % to ~ 4 % while for BDs it drops from ~ 60 % 
to ~ 33 %. This means the disc fraction drops by a factor 
~ 9 for LMS and by a factor of 2 for BDs in the period 
of time from ~ 3 Myr to ~ 7 Myr. This strongly suggest 
that, regardless of the initial disc fractions, disc evolution 
occurred faster for LMS than for BDs, which leads us to the 
third and last scenario. 

(Hi) The BD and LMS populations are coeval and the 
time scale for disc evolution depends on the mass of the cen¬ 
tral object, being slower for BDs than for LMS. Our results 
support this scenario, wh ich has been previously suggeste d 
by fe.g.j iRiaz et al.l ll2012ll and lLuhman fc Mamaie^ (l2012l l. 

More spectroscopic and photometric optical and near- 
IR observations for LMS and BDs populations in star form¬ 
ing regions with ages around ~ 10 Myr old and older are 
needed in order to clearly constrain how the evolution of the 
discs proceeds and what is its dependency with the mass of 
the central object. In order to understand the physical sce¬ 
narios involved, modeling and simulations are also needed, 
as well as a consistent statistical estimation of the number 
fractions and their comparisons for different regions. 
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